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Abstract 

A variety of studies have modeled the physics of material deformation and damage as examples of generalized 
phase transitions, involving either critical phenomena or spinodal nucleation. Here we study a model for 
frictional sliding with long range interactions and recurrent damage that is parameterized by a process 
of damage and partial healing during sliding. We introduce a failure threshold weakening parameter into 
the cellular-automaton slider-block model which allows blocks to fail at a reduced failure threshold for all 
subsequent failures during an event. We show that a critical point is reached beyond which the probability 
of a system-wide event scales with this weakening parameter. We provide a mapping to the percolation 
transition, and show that the values of the scaling exponents approach the values for mean-field percolation 
(spinodal nucleation) as lattice size L is increased for fixed R. We also examine the effect of the weakening 
parameter on the frequency-magnitude scaling relationship and the ergodic behavior of the model. 
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1. Introduction 

The physics of progressive material damage and healing under load has been studied using models that 
incorporate both failure thresholds and load sharing. Models of this class include slider-block models and 
fiber-bundle models. In both models, external loads are applied to a "loader plate" to which the blocks or 
fibers are attached. 

One type of avalanche model, using sliding blocks have been used to study the physics of frictional 
sliding, earthquakes [TJ [2J [3J, neural networks [3], as well as other driven threshold systems in nature. 
A variety of phenomena have been studied using these models, including aspects of scaling and critical 
phenomena [6l [7J [9] , as well as nucleation near the spinodal. A primary focus has been placed upon 
understanding the origins of the frequency-size relation seen in earthquake systems [Tl [21 El HI [HI IH1 El IH1 IS] - 
Slider-block models are closely related to other types of models for distributed failure, including fiber-bundle 
models [101 Ell EH E3] and fuse models [H| for the study of material damage and degradation and forest 
fire models [15] , 

Models for damage usually involve a system with a large number N of brittle elements such as fiber- 
bundles [TTJ EH1 EZ]- Often the elements have a statistical distribution of strengths. An external load is 
applied to the system, and several failure modes are typically observed. In some models, the fibers are 
completely brittle, so that failure occurs when the load on the fiber reaches the fiber strength. In other 
models, the strength of the element can degrade and weaken with time when the element load reaches a 
weakening value, with the time to failure dependent on the amount by which the stress exceeds the weakening 
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value [H3 [TT] [T21 E^]. In these models, an upper brittle strength value allows immediate failure when the 
stress on the element equals or exceeds the brittle strength. In these models, a damage variable a can be 
defined as the fraction of elements or fibers that have failed. 

Load sharing is an important property of both damage models and slider-block models. In the classical 
hierarchical fiber-bundle models [OH EHl EI] > all fibers are connected to end plates to which the external 
tensional stress is applied. When a fiber fails, its load is transferred in equal parts to all the other intact 
fibers through the end plates. Other details of fiber-bundle models have been discussed extensively in the 
literature (22J EH EH E3 E3 EH EH] . 

In the slider-block model, a loader plate is connected by springs to an array of blocks sliding on a 
frictional surface, with the blocks interconnected by coupling springs [HE]- The loader plate is typically 
driven either by a constant applied force F or by a constant applied velocity V. 

In the earliest slider-block models, a failure threshold (static friction threshold) was prescribed along 
with a residual or arrest stress level. Stress transfer in a slider-block model can occur through the coupling 
springs as well as through the loader plate. 

Healing is another physical process that occurs in some models of recurrent damage [29] . and can be 
observed in the laboratory and in nature. When an element fails, it transitions suddenly to a new state 
("ruptured state"). For a process such as cracking under tensile load, separation of the crack surfaces 
generally precludes healing. However, for shear failure, the crack surfaces typically remain in close proximity, 
and healing re-strengthening is possible under certain load conditions [29]. This is the case, for example, 
in sliding friction experiments, where strength increases as log(Ai) following a slip event, where At is the 
time interval since the slip. Presumably, the same process operates on natural faults and shearing surfaces, 
since sudden stick slip events are observed to reoccur regularly on these surfaces as earthquakes [301 [HI] , 
and therefore healing must occur. 

To summarize our results: We study a modified slider-block model for damage and failure in materials 
leading to a family of scaling exponents that can be measured in simulations. Our model differs from a 
traditional slider-block model in that we have included failure threshold weakening. Subsequent failures of a 
block during an event, will occur at a reduced failure threshold, parameterized by w. The failure threshold 
weakening parameter w represents the percentage reduction of the failure threshold after a block's initial 
slip. We investigate the behavior of the model as w is varied and show that above a critical value, w c , 
system-wide events occur regularly. The probability of a system-wide event scales with the distance from 
the critical point, w — w c . We examine the finite size effects and show in the limit of L — ¥ oo the scaling 
exponent /3, which characterizes the probability of system- wide events, approaches the corresponding order 
parameter exponent of mean-field percolation. 

This paper is organized as follows. In section 2 we discuss details of the model we use for simulations. In 
section 3 we present the results of our simulations, and in section 4 we discuss the mapping to percolation 
theory. In particular, we compare values of scaling exponents we obtain from the simulations to values for 
mean-field percolation and spinodal nucleation. We also place our results into the context of other recent 
work on similar models, with particular attention to the question of punctuated ergodicity [321 [SI 1331 134) 
which has been observed in both models as well as in natural earthquake fault systems. 

2. Model Definition 

The model we investigate here is a cellular automata slider-block model. We incorporate an idea for 
the physics of weakening and healing proposed in a different context |35j . which was a model for damage 
evolution in a continuum model of the earth's crust. Inasmuch as cellular automata slider-block models are 
simple systems in which to examine the physics of correlation and scaling, we have constructed a modified 
form of the usual slider-block model which has loader or pulling springs, connecting or interacting springs, 
and fixed failure and residual strengths. 

Our model includes loader and coupling springs as well as additional "weakening" and "healing" proper- 
ties |35j . By suitably tuning the threshold weakening, we find intervals in which a quasi-periodic earthquake 
cycle is observed that are dominated by large "characteristic" or "nucleation" events interspersed with time 



2 



Figure 1: (a) A sequence of event areas A are given as a function of event number (natural time) for zero weakening (ui = 0). 
A small segment of the total sequence is shown. The largest events are semi-periodic but are much smaller than the lattice area 
L 2 = 512 2 = 262, 144. The other parameters of the model run are: a F = 100, cr R = 25±25, R = 15 (q = 960), K C /K L = 0.10. 
(b) Frequency-magnitude statistics for weakening parameter w = 0. The frequency of events of area A, is plotted as a function 
of area A. The smaller events correlate well with the mean-field scaling relation taking r — 1 = 1.49 ± .05. 



intervals in which a Gutenberg-Richter, or scaling distribution of events is observed. In the following, we 
investigate and extend this model further, focusing in particular on the values of scaling exponents that are 
observed. It will be seen that the amount of threshold weakening can be described by means of a parameter 
w which appears to act as a relevant scaling field in the sense of nucleation and critical phenomena. 

Many dynamical models, for example forest fires |15j . display clusters of sites that are typically mapped 
onto percolation clusters, therefore our focus in this paper is on determining how the scaling exponents for 
the slider-block clusters are related to the scaling exponents typical of percolation theory. In particular, as 
we are interested in slider-block models with long range interactions, approaching mean-field, we examine 
mean-field percolation models in spatial dimensions d > 6 [36] . the upper critical dimension. 

An example of mean-field percolation is provided by the Bethe lattice (Cayley Tree), which is equivalent 
to percolation in d — > oo [36]). In that problem, two of the critical exponents have the values (Stauffcr 
and Aharony [see 36, table 2]): Order parameter exponent (3 = 1, and Fisher exponent r = 2.5. The 
remaining critical exponents can be recovered from the well-known scaling relations [36 . We note that for 
site percolation, the order parameter P(p) can be defined as the probability that a spanning or percolating 
cluster is present, i.e., a cluster whose bonds connect all sides of the lattice. 

Scaling in the slider-block model: In the model we consider here, each of the individual avalanches or 
clusters of "microscopic" blocks, which can be quite large, is considered to be a cooperative "macroscopic" 
slip event, or "earthquake" . Many studies have shown that if the frequency of events is plotted as a function 
of the number of event size (e.g. blocks in a cluster), the result is a power law of the form [5J [T51 |3"T1 155] : 

N=-^ T e-* Aa (1) 

where r and a are critical exponents. Here N is the number of clusters in the simulation with area A and 
e is a scaling field. In percolation theory, equation is called the Fisher- St auffer droplet model [36], and 
e oc (p — p c ) , where p is the occupation probability and p c is the value at the percolation threshold (critical 
point). In earthquake seismology, a similar scaling equation is called the Gutenberg-Richter frequency- 
magnitude relation [T] [21 SI El IFJ ■ The dynamics of the model are therefore reflected in the clusters of failed 
blocks and their statistics. 

2.1. Base Model 

We first introduce the standard cellular automata slider-block model. Our model consists of a 2d lattice 
of blocks sliding on a frictional surface. Our simulations are carried out on LxL square lattices ranging in 
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size from L = 256 to L — 1024, using periodic boundary conditions. Each of the blocks is coupled to a 
loader plate by means of a loader spring having spring constant Ki,. Each block is connected to q neighbor 
blocks by means of coupling springs having uniform spring constants Kc- The "total spring constant" is 
qKc + K L . The neighbor blocks form a square region with sides of length 2R + 1. Thus q = (2R + l) 2 — 1, 
where R is the range of interaction. Under rather general assumptions, it can be shown that the system can 
be regarded as approaching long range interactions R — > oo, the mean- field limit, as qKc — > oo [57]. Typical 
slider-block model dynamics occur when the loader plate is advanced in displacement until the force or stress 
(j on a single block meets or exceeds a failure threshold value ap, at which point the block is advanced a 
slip distance corresponding to a decrease in stress to the residual value. The residual stress value is drawn 
from a uniform distribution centered at or with a spread of San f° r each slip event. This small random 
component added to the slip distance has the effect of thermalizing the system [3^1 ■ 

Once the first block slips, stress is transferred to its q neighbors by means of the coupling springs, and 
in turn, other blocks may slip as well. The result is a cascade of slipping or failing blocks, often called 
an avalanche or event. We note that we use "zero- velocity" dynamics, meaning that the loader plate is 
advanced just enough so that only a single block is at the failure threshold at the initiation of an event and 
the loader plate is held fixed in location until the event has ceased. The loader plate is then advanced until 
the next event begins and the process is repeated. These events share much in common with avalanche 
phenomena in other physical systems, including sandpiles, neural networks, and electronic devices such as 
diodes. 

We illustrate the behavior of our base model in Fig. [lj In Fig. [lja) we plot the size of sequential 
avalanches A as a function of event number (natural time) . The size of an avalanche is characterized by the 
total area (number of blocks) that slips during the event. We give results for 10, 000 events from the middle 
of a simulation that included 10 6 events. In this simulation the system size is L = 512, the interaction range 
is R = 15 (q = 960), Kc/Kp = 0.1, ap = 100 and <jr — 25 ± 25. We see a wide range of rupture areas, 
the largest events occurring at relatively regular intervals. No system-wide events occur and the largest 
events, A w 25000 are much smaller than the system size L 2 — 262, 144. In Fig.JlJb) we plot the probability 
density function for the distribution of event sizes over the whole simulation run, the Gutenburg-Richter 
distribution. The frequency of events of area A is plotted as a function of area A. By suitably tuning 
the system parameters Kc/Kl and R the model shows a mean-field scaling relationship with exponent 
r- 1 = 1.5 j37j. 

2.2. Threshold damage model 

Our damage model uses the cellular automaton dynamics with a failure threshold ap for the initial 
failure and slip of all blocks and a residual stress an drawn from a random uniform distribution, as described 
above. We include damage into our model by considering any block that fails during an event to be weak or 
"damaged" for the remainder of the event. The first failure of all blocks occurs at an initial failure threshold 
of o~f. Any subsequent failure of a damaged block (one which has slipped once or more during the current 
event) will occur at a reduced failure threshold value ap * (1 — w), where < w < 1. The amount of damage 
to a block is controlled by the input weakening parameter w. Following termination of the avalanche of 
slipping blocks, the initial failure stress threshold of each block is reset to the original value ap. Thus we 
have damage, or weakening of the system following initial slip of each block by the ratio w, followed by 
healing back to the original failure threshold a F at the termination of the avalanche (macroscopic sliding 
event). In the following results, we investigate how the weakening parameter w effects the behavior of the 
system, particularly how beyond the critical value w c , the weakening parameter acts as a scaling field for 
the probability of system-wide events to occur. 

We next illustrate the behavior of our damage model when the weakening parameter w is greater than 
the critical value (w > w c ) so that system-wide events occur regularly. In Fig. [2](a) the rupture areas A are 
given as a function of event number (natural time) in Fig. [2] for w = 0.04 {w c = 0.012). System- wide events 
occur quasi-periodically in natural time. 

In Fig. [2](b) we give frequency-magnitude (Gutenburg-Richter) distribution for cluster sizes A for two 
values of the weakening parameter, w = 0.012 = w c and w — 0.04 > w c . The frequency of events of area A, 
is given as a function of area. For w — w c , shown with boxes, there are no system-wide events. For w > w c , 
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Figure 2: (a) A sequence of event areas A are given as a function of event number N (natural time) for weakening above 
the critical value (w > w c ). A small segment of the total sequence is shown. System-wide events A = L 2 = 512 2 = 262, 144 
occur regularly. The other parameters of the model run are: ap = 100, <tjj = 25 ± 25, R = 15 (q = 960), Kq/Kl = 0.10. 
(b) Frequency-magnitude statistics for two values of the weakening parameter: ui = w c and w > w c . The frequency of events of 
area A, is plotted as a function of area A. At the critical point ui = w c = 0.012 there are no system-wide events and the smaller 
events correlate well with equation[l]taking r — 1 = 1.47 ± .05. Above the critical point w = 0.04 there are system-wide events 
occurring regularly producing the spike at A = 262, 144, while the smaller events continue to obey the GR scaling relationship 
with r very near 2.5. 



shown with circles, system-wide events occur regularly, producing the large spike at area A = 262, 144. The 
system-wide events have a different structure than the smaller events. The smaller events have a diffuse 
structure with a mixture of failed and unfailed blocks where each failed block typically slips only once during 
the event. The system- wide events are dense with all blocks failing and each failed block is likely to slip 
multiple times during the event. Also included in this figure is the correlation between the w = w c data 
and the straight-line relation given in equation [l] Good agreement is obtained by fitting the region from 
A = 1 — 1000 giving a scaling exponent r — 1 = 1.47 ± .05. The system parameters are the same as in[T] 
Similar results are obtained for simulations using a larger lattice sizes L = 768 and 1024. 



3. Simulation Results 

The primary purpose of our simulations is to study the role of the weakening parameter w. All simulations 
will assume a failure stress op — 100, a residual stress or = 25 ± 25, and a loader spring constant Kl = 1 
and periodic boundary conditions. For no weakening (base model), w = 0, system-wide clusters of failed 
sites in slider-block models are exceedingly rare, and appear only in the mean-field limit as qKc °°i 01 
alternatively as the range of interaction R — > oo |40) . In our simulations, we find that as the weakening 
parameter w increases (w > 0), a value is reached at which system- wide events ("characteristic earthquakes") 
begin to appear. Indeed, this result was noted for a different model in 35 , which demonstrated that it is 
possible to tune such a model so that intermittency in the appearance of system-wide avalanches is observed. 
This suggests that we view the existence of system- wide events of failed sites in our slider-block model in a 
similar way to the infinite cluster in percolation, namely, that system-wide events should be associated with 
an order parameter and that there are associated critical exponents. 

Since the weakening parameter w controls the nature of the scaling in this model, we conjecture that w 
plays a similar role to the occupation probability p in site percolation. Thus we define an order parameter 
P(w) to be the ensemble average of the probability of a site (block) selected at random being part of a 
system- wide cluster ( "characteristic event" ) . Here we consider system- wide events to be those that have a 
total area equal to the lattice area. The order parameter we have defined is analogous to the 'strength of 
the infinite network' order parameter P(p) defined in percolation theory |36) . By defining a "characteristic 
events" as encompassing the entire lattice area, we can calculate P(w) simply as the fraction of observed 
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Figure 3: The order parameter vs the weakening value adjusted by the critical weakening value. The order parameter is defined 
as the probability of of a system wide event occurring. The plot is on a logarithmic scale emphasizing the power-law dependence 
of P(w — w c ) with a scaling exponent very near = 1. The parameters of the model are the same as in Fig. [T] 

avalanche events that are system-wide events during the simulation run. In general we find that for w < w c , 
P(w) = in any finite time interval, but that for some value, P(w) > 0. 

The value of w c is determined numerically and is equal to the value of w beyond which system- wide events 
("characteristic earthquakes") begin to appear. We have determined the critical value of the weakening 
parameter w c when system-wide events ("characteristic earthquakes") begin to appear is w c = 0.012 for the 
simulation parameters used in Figs. [I] and [2] As discussed above we define the order parameter P(w — w c ) 
as the ensemble average of the probability a site selected at random is part of the system- wide event, where 
a system-wide event encompasses then entire lattice. We calculate the order parameter as the fraction of 
observed avalanche events that are system-wide events. Near a critical point we would expect to find that: 

P(w - w c ) oc(w- w c f. (2) 

In Fig. [3] we plot the dependence of logP on log(u> — w c ) for the same system parameters of the model 
used in Fig. [T] Each data point represents a simulation run with differing value of w. The range of w is 
bounded on the right by the system that has all events spanning the entire lattice, that is P(w — w c ) = 1. 
The value of scaling exponent for the initial scaling region is near 1, the same value that characterizes mean- 
field percolation as defined in equation [2] The observed j3 = 0.94. As we increase the system size, smaller 
values of w yield non-zero probabilities of system-wide events, extending the scaling region to the left. In 
addition, the larger system sizes produce a slightly larger value of the scaling exponent /3 with an apparent 
asymptotic value of /? = 1. We have measured the scaling exponent for systems with lattice dimension of 
256,512 and 768. 

4. Analysis 

4-1. Scaling in non- equilibrium systems 

It is well known that high dimensional models composed of interacting sites, spins, trees, cells, units, 
or blocks demonstrate clustering and scaling phenomena [HI 02j [43J S3]- Ising models, which are often 
used to demonstrate both equilibrium phase transitions (second order) and non-equilibrium transitions 
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(nucleation), have been extensively analyzed |41j . More recently, it has been observed that non-equilibrium 
systems, including sandpile (self-organized critical) and forest fire models [43j [44] also exhibit clustering 
and scaling phenomena. For the case of slider-blocks, it was found that the energetics of the model in the 
near mean- field regime lead to a Boltzmann factor for the energy fluctuations [39, 8, 33 . The Boltzmann 
distribution of energy fluctuations, along with the apparent punctuated ergodicity as measured by the TM 
metric, described below, imply that the scaling properties of these models may have a similar origin to 
scaling properties in equilibrium systems. Later, it was found that the same ergodic property holds for 
coupled map lattices in the near mean- field regime |45j . 

For Ising models, it has been known for many years [46] that the scaling exponents associated with 
the Ising critical point can be mapped onto a percolation model constructed from a correlated site-bond 
probability. By themselves, random-site percolation clusters are not isomorphic to Ising clusters, because 
Ising clusters have both a correlated component, as well as a random (geometric) component. However, 
it was shown |46) that a bond probability can be constructed to define bond percolation clusters from the 
Ising spin clusters that are in fact isomorphic to random site percolation clusters. The critical exponents are 
the same, and the percolation critical point maps onto the Curie temperature defining the ferromagnetic- 
paramagnetic transition. Thus the bond probability is used to separate the correlated and geometric effects 
that define the Ising clusters. 

For slider-block models, the physics is different but related, and in fact, simpler. Here a block fails only if 
the stress is large enough. So as Klein et al. [37l pg. 64] note, all clusters of failed blocks are constructed from 
correlated components, and there is no need for an additional bond probability. For that reason, we expect 
that the scaling exponents should map directly onto those for an associated site percolation problem. Since 
we are dealing with a near mean-field model having long range interactions, this implies that the scaling 
exponents for the slider-block model should be the same as the Bethe lattice (mean-field percolation), as the 
range of interaction of the slider-block model becomes large. Determining the validity of this correspondence 
is a major motivation for the results presented here. 

We have shown in Figs.Qb) and[2jb), the frequency-magnitude scaling for the slider-block model obeys 
the scaling relation given in [T] In Fig. [3] we have also indicated the scaling relation between the order 
parameter P and the scaling field w — w c introduced into our damage model. The values r — I = 1.51 and 
/3 = .94, provided we are using a finite lattice, are within experimental error of the corresponding mean-field 
percolation model scaling exponents. 

4-2. Ergodic Dynamics 

The slider-block model is an example of a driven nonlinear threshold system where the dynamics of 
the system are strongly correlated in space and time. It has been shown that as the range of interaction 
becomes large, the slider-block model approaches a mean-field spinodal resulting in scaling phenomena, such 
as the frequency-magnitude (Gutenburg-Richter) shown in Figjlja) [37J. In this limit, the system resides in 
a metastable equilibrium, where equilibrium-like properties, such as stationary dynamics, equally probable 
microstates, and effective ergodicity all hold true. This metastable state is disrupted by large events, and 
the system evolves to a new metastable state. Therefore, the slider-block model displays punctuated or 
intermittent ergodicity [4"T] . 

A measure of ergodicity can be obtained using the Thirumalai-Mountain (TM) fluctuation metric [3"21r4"5] . 
The TM metric, f2(i), measures the difference between the time average of an observable at each site and its 
ensemble average over the entire system. The TM metric relies on the idea of statistical symmetry, where 
the statistics of one particle (block) is presumed identical to those of the entire system. Here the observable 
is the stress on each block a(t). The metric is defined as follows: 

i 

where N is the number of blocks, and the quantities (Ji{t) and (cr(t)) are given by 

<7i(t) = i f dt'atit') (4) 
{ Jo 
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Figure 4: The inverse TM metric for measuring effective ergodicity is overlayed with the event size sequence as a function of 
loader plate position, for the "threshold damage model" at the critical point w = w c . A system is effectively ergodic if the 
stress metric decreases as l/t. In (a) we plot the inverse TM metric vs loader plate position (time). This plot shows a linear 
trend, displaying the system is effectively ergodic from Vt = 300 — 600. In (b), we show a detailed view of a break in the 
ergodicity caused by a large event. 

and 

(*(*)> = ^ !>(*)■ ( 5 ) 

i 

If the system is effectively ergodic at long times, the TM metric decreases in time as, fl(t) = ^ [35], where 
D e is a diffusion constant related to the rate at which the phase space is explored. In Fig. [4] we plot the 
inverse of stress metric introduced in [37]. As noted above, a system is effectively ergodic if the stress metric 
decreases as l/t, therefore, a plot of the inverse metric, ^, vs time, t, will show a linear relationship for 
ergodic systems. In Fig. Qa), we plot the inverse stress metric, with "boxes", and the sequence of events, 
with "dots", as a function of the loader plate position (time) for the "threshold damage model" at the 
critical point (w — w c ). The inverse TM metric shows a linear trend from loader plate position 300 — 600. 
The initial curved portion of the plot is due to transient effects. In Fig.[4|b), we show a small section of the 
plot in Fig. Qa) centered on a large event of size A ~ 50000. Here we show that large events disrupt the 
ergodicity of the model for a short period of time as evident by the "kink" in the inverse TM metric. This 
model has punctuated ergodicity. In Fig. [5] we plot the inverse metric for the "threshold damage model" 
above the critical point. Here the system is not ergodic over any time period. The grid size events prevent 
the system from residing in a metastable equilibrium for a long period of time. The inverse metric is well fit 
to a quadratic function as noted in |32j for systems that are non-ergodic. This result is in agreement with 
a different adaptation of "damage" into the long-range slider-block model proposed by [34]. Here blocks in 
the lattice are removed and no longer contribute to the dynamics of the system after a specified number of 
failures. At the onset of the blocks becoming "damaged" the system loses its ergodicity. 

5. Discussion 

5.1. Application to Data 

It is useful to compare the scaling exponents determined from our numerical simulations with scaling 
exponents seen in laboratory experiments and in nature. In all of these experiments, the fracturing and 
acoustic emissions are mixed-mode, not simple shear sliding as is the case for our simulations. As a result, we 
expect to observe a range of values for scaling exponents. Typically, the data are in the form of probability 
density function (PDF) as a function of energy E of acoustic emissions (AE). For frequency vs. energy 
measurements, and assuming that area of microcrack is proportional to energy released, which would be the 
case if crack area is proportional to crack energy for these microcracks, a range of scaling exponents r — 1 
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Figure 5: The inverse TM metric and the event size sequence are shown for the "threshold damage model" above the critical 
point (w > w c ). The inverse metric shows a 1/t 2 trend, displaying that this model is not ergodic over the observed time period. 
The grid-size events, shown with "dots", prevent the system from residing in a meta-stablc equilibrium for any extended period 
of time. 

are observed, from 1.2 to 2.0, with most values clustering around 1.5 as is observed in our simulations 051 
[SOI [STJ [55] . On the other hand, if crack energy released is proportional to crack area to the 3/2 power, as is 
the case for macroscopic cracks, then we would expect scaling exponents at the lower end of the observed 
range, clustering around r — 1 = 1.0. 

5.2. Summary 

In this paper, we have discussed the various types of simulations used to model the physics of material 
deformation and damage, including fiber-bundle models and slider-block models. We have also discussed 
the experimental and observational motivation for incorporating damage and healing into these models. We 
infer that healing must occur on faults where stick-slip events (earthquakes) occur regularly, motivating the 
addition of damage and partial healing into a traditional earthquake model. We included the concept of 
damage and healing into the standard cellular-automata slider-block model by means of failure threshold 
weakening. We reduced ( "damaged" ) each block that slips during an event by allowing subsequent failures of 
the damaged blocks to occur at a lower failure threshold parameterized by w. After the event has terminated 
we heal all blocks to their original failure threshold levels. 

We have shown that the weakening parameter w behaves as a scaling field for the ensemble average 
probability of a site selected at random to be part of a system-wide event P(w — w c ), just as the occupation 
probability is a scaling field for the strength of the infinite network in percolation theory. We also discussed 
the mean-field percolation problem along with the various critical exponents that can be measured. The 
results of our simulation yield scaling exponents that approach those of the in Bethe lattice (mean-field) 
percolation problem, specifically (3 — 1. We have demonstrated that the "threshold weakening model" 
displays ergodic behavior up to the critical point. Beyond the critical point w c , the model is no longer 
ergodic. We conclude that "damage" as we have modeled it, is a non-equilibrium process. We have also 
shown the scaling exponents found in our slider-block model to be within the range of those found in 
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experiments. We conclude that the mean-field interactions may dominate the failure of material. Further 
studies of this model including an investigation of time dependent healing are currently being studied. 
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